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By applying a local Rashba spin-orbit interaction on an individual quantum dot of a four-terminal 
four-quantum-dot ring and introducing a finite bias between the longitudinal terminals, we theoret- 
ically investigate the charge and spin currents in the transverse terminals. It is found that when the 
quantum dot levels are separate from the chemical potentials of the transverse terminals, notable 
pure spin currents appear in the transverse terminals with the same amplitude and opposite polar- 
ization directions. Besides, the polarization directions of such pure spin currents can be inverted by 
altering structure parameters, i.e., the magnetic flux, the bias voltage, and the values of quantum 
dot levels with respect to the chemical potentials of the transverse terminals. 
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Since the original proposal of the spin field effect 
transistor by Datta and Das,[l[ enormous attention, 
from both experimental and theoretical physics com- 
munities, has been devoted to the controlling of the 
spin degree of freedom by means of the spin-orbit 
(SO) coupling in the field of spintronics. [2j Partic- 
ularly, in low-dimensional structures Rashba SO in- 
teraction comes into play by introducing an electric 
potential to destroy the symmetry of space inversion 
in an arbitrary spatial direction. [H, 0>S [H, 0] Thus, 
based on the properties of Rashba effect, electric 
control and manipulation of the spin state is feasible. 
Accordingly, the Rashba-related electronic proper- 
ties in mesoscopic systems have been the main con- 
cerns in spintronics, such as spin decoherencejl, [§] 
and spin current [H| • 

Quite recently, Rashba interaction has been intro- 
duced to coupled quantum dot (QD) systems. Be- 
cause the coupled QD systems possess more tunable 
parameters to manipulate the electronic transport 
behaviors, a number of interesting Rahsba-induced 
electron properties are reported, EH moreover, 
it is theoretically predicted that pure spin cur- 
rents are possibly realized in a triple-terminal QD 
structures only by the presence of a local Rashba 
interaction [TJ, Ell- Following such a hot topic, in 
this paper we propose a new theoretical scheme to 
realize the pure spin current by virtue of Rashba in- 
teraction. We introduce Rashba interaction to act 
locally on one component QD of a four-QD ring with 
four terminals. Our theoretical investigation indi- 
cates that the unpolarized charge current injected 
through the longitudinal terminals gives rise to the 
emergence of pure spin currents in transverse termi- 
nals with the same amplitude and opposite polariza- 
tion directions, and the polarization direction of the 



* Author to whom correspondence should be addressed. Email 
address: zys@mail.jlu.edu.cn 




-3-2-10123 

e (in units of r) 




>|> (in units of Jt) 



FIG. 1: (a) Schematic of a four-terminal four-QD ring 
structure with a local Rashba interaction on QD-2. Four 
QDs and the leads coupling to them are denoted as QD-j 
and lead-j with j = 1 — 4. The currents vs the QD levels 
Sj as well as the magnetic phase factor <f> are shown in (b) 
and (c), respectively. The parameter values are Fj = 2F 
and q = 0.4. In (b), Ej = eo and = and in (c), 
the QD level Sj — F. (d) The currents vs <j) in the case 



pure spin current in either terminal is tightly depen- 
dent on the adjustment of structure parameters. 

The structure under consideration is illustrated in 
FigOJa). In such a system four leads are coupled to 
their respective QDs in a four-QD ring, with the 
Rashba effect applied to act locally on QD-2. The 
single-electron Hamiltonian in this structure can be 
written as H s = -^+V r (r) + ^-[a(CTXp) + (<7Xp)a], 
where the potential V(r) confines the electron to 
form the structure geometry. The last term in H s 
denotes the local Rashba SO coupling on QD-2. For 



2 



the analysis of the electron properties, we select the 
basis set {il>jkX<r>' l l ) jXcr}(' , fij an d ipjk are t ne orbital 
eigenstates of the isolated QDs and leads in the ab- 
sence of Rashba interaction with j=l-4; Xa with 
a =1,1 denotes the eigenstates of Pauli spin oper- 
ator <x z ) to second-quantize the Hamiltonian, which 
is composed of three parts: H s = H c + H-d + Ht- 
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the transmission function, describing electron tun- 
neling ability between lead-j to lead-j', and fj(u>) 
is the Fermi distribution function in lead-j. Tj = 
2Tr\Vj a \ 2 pj(uj), the coupling strength between QD-j 
and lead-j, can be usually regarded as a constant. 
G r and G a , the retarded and advanced Green func- 
tions, obey the relationship [G r ] = [G a ]t. From 
the equation-of-motion method, the retarded Green 
function can be obtained in a matrix form, 
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where c^ ka and dj a (cjk. a and dj a ) are the creation 
(annihilation) operators corresponding to the ba- 
sis in lead-j and QD-j. Ejk and Sj are the single- 
particle levels. Vj a = (ipjXa\H\ipjkXa) denotes QD- 
lead coupling strength. The interdot hopping am- 
plitude, written as ti„ = t\ — iasi(l = 1,2), has 
two contributions: U = (ipi\HQ\tpi + i) is the ordi- 
nary transfer integral independent of the Rashba in- 
teraction; si = i(ipi\ap x + p x a\ipi+i), a real quan- 
tity for real ipi and indicates the strength of 
spin precession. Finally, in the Hamiltonian r; = 
(ipi\ap z +p z a\ipi+i) is a complex quantity represent- 
ing the strength of interdot spin flip. To get an in- 
tuitive impression about the typical values of these 
parameters in the Hamiltonian, we assume that each 
QD confines the electron by an isotropic harmonic 
potential |m*wor 2 - Then, the four QDs distribute 
on a circle equidistantly, and the interval in be- 
tween is 21q(Iq = y/h/m*LJo). Besides, we assume 
that the electron occupies the ground state in each 
QD. By defining a dimensionless Rashba coefficient 
as a = a / (Shuolo) , we obtain the rough relation 
of the parameters: t\ = t2, Si = S2, T\ = — T2, 
and \si\ = \r t \ ~ at t . Thereby we can express the 
interdot hopping amplitude in an alternative form: 
ti a = i/vl + a 2 e~ laip with ip = tan -1 a. Here the 
Rashba interaction brings about a spin-dependent 
phase factor, which can be tuned by varying the 
electric field strength. In the above Hamiltonian 
the phase factor attached to £4 accounts for the 
magnetic flux through the ring. In addition, the 
many-body effect can be readily incorporated into 
the above Hamiltonian by adding the Hubbard term 

V e -e = 'Ylja ~2 n j<y n jS- 

Starting from the second-quantized Hamiltonian, 
we can now formulate the electronic transport prop- 
erties. With the nonequilibrium Keldysh Green 
function technique, the current flow in lead-j can 
be written as [la] 
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In the above expression, gj a is the Green function of 
QD-j unperturbed by the other QDs and in the ab- 
sence of Rashba effect. gj a = \(z — ej)\j a 

with z = u + i0+ and \ ja = z _ £j .!^ j . +t ,. <nj ._ ) 
resulting from the second-order approximation of 
the Coulomb interaction 17 1. n,j a = (d-dj C ) can 
be numerically resolved by iteration technique with 
{n jtT ) = -£fduGf ajja , 

We now proceed on to calculate the currents in 
lead-j, lead-1 and lead-3 in this case. As for the 
chemical potentials in respective leads, we consider 
pi as the zero point of energy of this system and 
pi = /13; P2 and ^4, the chemical potentials in other 
two leads, are considered with P2 = pi + eV/2 and 
/i4 = p\ — eV/2, in which V is the bias voltage. The 
charge and spin currents are defined respectively as 

Jrnc = Jrn] + Jml an d J ms = J m | — J m | (m = 1, 3). 

Before calculation, we introduce a parameter T as 
the unit of energy, the order of which is meV for 
some experiments based on GaAs/GaAlAs QDs, as 
mentioned in the previous works \l§L [l9|. 

To carry out the numerical calculation about the 
spectra of the charge and spin currents in lead-1 and 
lead-3, we choose the Rashba coefficient a = 0.4 
which is available in the current experiment . pol| The 
structure parameters are taken as |t; CT | = t$ = = 
T, and the bias voltage is eV = 2T. In FigfJIb), 
the currents versus the QD levels are shown in the 
absence of magnetic field. It shows that in the case 
of the QD levels separate from the energy zero point 
(i.e., £o 7^ 0), besides the emergence of pure spin 
currents, a more interesting phenomenon is that the 
amplitude of J± s is almost the same as that of J3 S . 
However, in the case of £o > the value of J% s is 
greater than zero and Js s < 0, whereas J\ s < and 
J^ s > under the condition of £o < 0. This means 
that, in such a structure applying a finite bias be- 
tween the longitudinal terminals can efficiently bring 
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FIG. 2: The spectra of transmission functions 
T ma>m / cr (m=l,3 and m'=2,4). (a) and (b) Zero mag- 
netic field case, and (c)-(d) magnetic phase factor <j> = ^. 
The domains of integration to calculate the currents are 
labeled by the arrows. 



out transverse pure spin currents with the opposite 
polarization directions of them. Moreover, when the 
QD levels exceed the zero point of energy of the sys- 
tem the polarization directions of these pure spin 
currents will be thoroughly changed. However, for 
the case of the QD levels consistent with the zero 
point of energy, e = 0, no pure spin current comes 
about despite a nonzero magnetic flux through the 
QD ring, as shown in Figfljc). When the QD lev- 
els take a finite value (eg., £o = T), as shown in 
FigOJd), not only there are apparent spin currents 
in the transverse terminals (lead-1 and lead-3), but 
also with the adjustment of magnetic flux in cither 
transverse terminal the charge and spin currents os- 
cillate out of phase. In the vicinity of cf> — (n— 
J\ c and Js c reach their maxima; Simultaneously, the 
spin current J ms are just at a zero point. On the 
contrary, when <f> = nix the situation is just inverted, 
the maximum of J ms encounters the zero of J mc . In 
particular, with the change of magnetic phase fac- 
tor from (f> — 2nir to <fi — (2n + 1)tt the polarization 
directions of the transverse pure spin currents are in- 
verted. So, it should be noticed that tuning the QD 
levels to an nonzero value with respect to the zero 
point of energy is a key condition of the appearance 
of transverse pure spin currents. 

The calculated transmission functions are plotted 
in Fig. [2] with £j = T. They are just the integrands 
for the calculation of the charge and spin currents 
(see Eq.Q). By comparing the results shown in 
Figs|2ja) and[2h), we can readily see that in the 
absence of magnetic flux, the traces of Tit,2Ti ^i|,4i> 
2^31,24.5 and 2~3T,4T coincide with one another very 
well, so do the curves of Tij^j., ^1T,4T> ^3T,2T; and 
Tzi^i- Substituting such integrands into the cur- 
rent formulae, one can certainly arrive at the result 
of the distinct pure spin currents, which flows from 



lead-1 to lead-3 in such a case. On the other hand, 
these transmission functions depend nontrivially on 
the magnetic phase factor, as exhibited in FigJ^c) 
and (d) with <fi = ^. In comparison with the zero 
magnetic field case, herein the spectra of T mr j^„ are 
reversed about the axis lu = T without the change 
of their amplitudes, but T ma ^ a only present the en- 
hancement of their amplitudes. Similarly, with the 
help of Eq.Q, one can understand the disappear- 
ance of spin currents in such a case. 

The underlying physics being responsible for the 
spin dependence of the transmission functions is 
quantum interference, which manifests if we analyze 
the electron transmission process in the language 
of Feynman path. Notice that the spin flip arising 
from the Rashba interaction does not play a leading 
role in causing the appearance of spin and charge 



currents [15j. Therefore, to keep the argument 
simple, we drop the spin flip term for the analysis of 
quantum interference. With this method, we write 
Tia,2a = |tict,2ct | 2 where the transmission probability 
amplitude is defined as ti^o- = Vi a G r 1<7 2cr V2<T with 

Vj a = Vj a ^2Trpj(u>). By solving G r la2 a-> we nn< i 
that the transmission probability amplitude T\ a ^a 
can be divided into three terms, i.e., Tio- 2a = t 
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h^lcr9l<7 t 4 e ~ l ' l '94<T t 393at2a92aV2a, and 
~T5 V *a 9l<7tla92aha93at2a 92aV 2a with 

D = det{[G r ] _1 }rij5^- % observing the 
structures of TfJ M , t^ 2 „, and t^ 2 „, we can 
readily find that they just represent the three 
paths from lead-2 to lead-1 via the QD ring. 
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' ia.2a 1S ^^20- = [(/>—2a(p\. It is clear that only these 
two phase differences are related to the spin polar- 
ization. Tict.40- can be analyzed in a similar way. We 

then write T\ a ^ a 
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difference between t\ 1 J a „ and 
2cnp + 
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and 
The 
„ is 



[<f>-2a V ] 



A^J = [<f>- 2a<p + 2 + 3 ] , and A^fj 
originates from the phase difference between r}^ 4(7 
and t^ Acf . Utilizing the parameter vlues in Figj2l 
we evaluate that <p « j and 9j = — ^ at the point of 
lu = 0. It is apparent that when = only the phase 
differences Atfr^J and A</> 4 ^ are spin-dependent. 



Accordingly, we obtain that At/)^ 



and A^IV = —AifiVJ = — 4? , which clearly prove 
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that the quantum interference between t{^ 2 | and 
T iT^2T ( r U^4i ano " T U 41 alike) is destructive, but the 



4 



constructive quantum interference occurs between 

r^J 21 an d 2J. ( T it 2T an< ^ T 1T 2T a like) . Then 
such a quantum interference pattern can explain 
the traces of the transmission functions shown in 
Fig[2ja). In the case of <p = \^ we find that only 

(2) 

A^ga),, are crucial for the occurrence of spin polar- 
ization. By a calculation, we obtain A^S^j = ff 
and A^S^j = ff > which are able to help us clarify 
the results in Figl^c) and (d). Up to now, the 
characteristics of the transmission functions, as 
shown in Fig(2l hence, the tunability of charge 
and spin currents have been clearly explained by 
analyzing the quantum interference between the 
transmission paths. 
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FIG. 3: The currents versus eo in the presence of many- 
body terms with Uj = 2F and 3r, respectively. 

So far we have not discuss the effect of electron 
interaction on the occurrence of pure spin currents, 
though it is included in our theoretical treatment. 
Now we incorporate the electron interaction into the 



calculation, and we deal with the many-body terms 
by employing the second-order approximation, since 
we are not interested in the electron correlation 
here. Figure [3] shows the calculated currents spectra 
with Uj = 2T and 3r, respectively. Clearly, within 
such an approximation the Rashba-induced trans- 
verse pure spin currents remain, though the current 
spectra oscillate to a great extent with the shift of 
QD levels. 

In conclusion, by introducing a local Rashba in- 
teraction on an individual QD, we have studied the 
electronic transport through a four-QD ring with 
four terminals. As a consequence, the Rashba- 
induced transverse pure spin currents are observed 
by applying a finite bias on the longitudinal termi- 
nals. The modulation of the QD levels and the mag- 
netic phase factor can efficiently adjust the phases 
of the transmission paths, thus the spin-dependent 
electron transmission probabilities can be controlled 
by tuning these structure parameters, which brings 
out the change of the amplitudes and directions of 
the pure spin currents. With respect to the quan- 
tum interference in such a structure, we have to il- 
lustrate two aspects. First, the applying of Rashba 
interaction is the precondition of the spin-dependent 
electron transmission. The presence of the multi- 
terminal configuration leads to the comparative am- 
plitudes of the different transmission paths for the 
quantum interference. Finally, it should be empha- 
sized that altering the longitudinal bias, equivalent 
to interchange the sequence numbers of lead-1 and 
lead-3, can also change the polarization directions of 
the pure spin currents. 
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